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Observational effects from quantum cosmology 
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The status of quantum cosmologies as testable mod- 
els of the early universe is assessed in the context of 
inflation. While traditional Wheeler-DeWitt quantiza- 
tion is unable to produce sizable effects in the cosmic 
microwave background, the more recent loop quan- 
tum cosmology can generate potentially detectable 
departures from the standard cosmic spectrum. Thus, 
present observations constrain the parameter space 
of the model, which could be made falsif iable by near- 
future experiments. 
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During the last years, quantum gravity has been receiv- 
ing a great amount of attention from the community 
of theoreticians. The driving motivation, familiar to any- 
one who has tried his or her fortune at least once in 
this broad subject, is to realize a consistent, ultraviolet fi- 
nite merging of general relativity with quantum mechan- 
ics. The programme can be carried out in various forms, 
from ambitious theories of everything (such as string the- 
ory) where all forces are unified to more minimalistic 
approaches aiming to quantize gravity alone. In the lat- 
ter category fall loop quantum gravity (TQG), asymptotic 
safety, spin-foams, causal dynamical triangulations and 
many others [1]. 

A problem endemic to most of these scenarios is 
their difficulty in making contact with observations. This 
stems from the highly technical nature of the theoreti- 
cal frameworks, where the notions of conventional geom- 
etry and matter, continuum spacetime, general covari- 
ance and physical observables are typically deformed, 
modified, or disappear altogether. The lack of exper- 
imental feedback makes it quite difficult to discrimi- 
nate among different models and, chiefly, to characterize 
them as falsifiable. 

It is natural to turn to cosmology in an attempt 
to bridge this gap and advance our knowledge [2, 3]. 
The early Universe is an ideal laboratory where extreme 



regimes of high energy and high curvature are realized. 
Under such conditions, it is expected that quantum grav- 
itational effects become sizable. Also, the symmetry re- 
duction entailed in cosmological settings decimates the 
degrees of freedom of background-independent theories 
and allows one to simplify the latter to a technically man- 
ageable level. The resulting models retain some (or most) 
of the main features of the full theory and can be better 
manipulated to extract observables. 

Canonical quantum gravity is a popular example of 
this mechanism. The present review focusses on two of 
its incarnations, namely, the traditional Wheeler-DeWitt 
(WDW) model (e.g., [4,5]) and the more recent loop quan- 
tization [6] . The most ancient phase about which we have 
gathered experimental data is inflation, a period of accel- 
erated expansion of the universe which left a relic in the 
cosmic microwave background (CMB) radiation. A study 
of the inflationary perturbations and the associated spec- 
tra allows us to track down quantum corrections and con- 
front them with the observed CMB power spectrum. Al- 
though the outcome of this procedure is a constraint on 
the free parameters of the models rather than an actual 
prediction, time seems ripe for the very next generation 
of experiments to exclude notable portions of parameter 
space. As a minimal present-day achievement, we can at 
least state that quantum cosmology models are compati- 
ble with observations. 

The stark contrast between the type of quantum cor- 
rections arising in these scenarios highlights how sensi- 
tive the physics is of the quantization scheme and vari- 
ables. The typical energy scale during inflation is es- 
timated to be about the grand-unification scale, H ~ 
10 15 GeV, corresponding to an energy density p in fl ~ 
H 2 ie\ x ~ 10 68 GeV 4 .HereH:= a I a is the Hubble parame- 
ter, a is the scale factor of the universe and a dot denotes 
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differentiation with respect to synchronous time. In con- 
trast, classical gravity is believed to break down at dis- 
tances shorter than the Planck length £p\ = V Gh, i.e., at 
energies above 10 19 GeV. The ratio between the inflation- 
ary and Planck energy density is very small, 



^ ~ {e Pl H) 2 ~ 10- 

Ppi 



(1) 



and quantum corrections are expected to be of the same 
order of magnitude or lower. Thus, quantum-gravity ef- 
fects would be, in fact, well below any reasonable experi- 
mental sensitivity threshold, at least as far as inflation is 
concerned. WDW quantum cosmology realizes precisely 
this type of corrections and endorses the above naive ar- 
gument. 

On the other hand, the polymeric quantization of 
loop quantum cosmology (LQC) [5,7,8] generates correc- 
tions which are not of the form (1). To get a rough idea 
of how these corrections arise, one begins by observing 
that geometry operators representing areas and volumes 
acquire a discrete spectrum in this context. This is be- 
cause states of loop quantum gravity, spin networks, are 
graphs whose edges e are labeled by quantum numbers 
j e . An area intersected by some of these edges is deter- 
mined by these quantum numbers, giving the spectrum 
sd = y^Ze \j jeije + 1), where y < 0(1) is the Barbero- 
Immirzi parameter. One single edge defines an "elemen- 
tary plaquette" of area oc (^y/ j e (je+ 1); the latter fea- 
tures the Planck area but its actual value depends on 
the spin quantum number. Since calculations on realistic 
graphs are very hard in the full theory, it is convenient to 
focus one's attention on a simplified phenomenological 
setting. In particular, a homogeneous quantum inflation- 
ary universe with small inhomogeneous perturbations 
may be represented by a quantum semi-classical state 
characterized by a length scale L. This scale is thought 
of as encoding the discreteness of the geometry. Any re- 
gion of volume T = a 3 Y (arbitrary, if spatial slices are 
non-compact) can be decomposes into discrete patches 
of size ~ L 3 . The inflationary scale is thus replaced by an 
effective quantum-gravity scale 



(2) 



In general, inverse powers of L cannot be quantized to 
a densely defined operator because the spectrum of the 
volume contains 0. This requires to reexpress inverse vol- 
umes (appearing in the constraints and hence in the 
dynamics) via Poisson brackets, which in turn feature 
derivatives by L. Quantum discreteness then replaces 
classical continuous derivatives by finite-difference quo- 
tients. For example, the expression {2\/~Ly l = d\f~LldL 



would become {\/L + £p\- \/L-£p\)/[2£pi), strongly dif- 
fering from {2\fTy l when L is as small as the Planck 
length, L~ lp\. For larger L, corrections are perturbative 
and of the order (p\IL, so in general the type of inverse- 
volume quantum corrections are expressed by the ratio 



Ppi \ L ) ~ 



(3) 



In practice, the actual size of LQC effects will lie well be- 
low the over- optimistic upper bound (3), but above the 
naive estimate (1). 

After introducing the theoretical frameworks in sec- 
tions 2.1 and 3.1, CMB observations will be used to 
pin down these effects (sections 2.2 and 3.2). For the 
Wheeler-DeWitt model, we shall do so in considerably 
more detail than can be found in the present litera- 
ture; section 2.2 contains original material. The scantly 
touched topic of non-Gaussianity in quantum cosmol- 
ogy will be also discussed (section 4). In the following, 
H = l = c. 



2 Wheeler-DeWitt cosmology and 
observations 

2.1 The model 

2.1.1 Homogeneous background 

In canonical formalism, symmetry and dynamics are en- 
conded in a set of constraint equations valid on dynam- 
ical trajectories. For gravity and matter, the total Dirac 
Hamiltonian [9] obtained after imposing second-class 
constraints and skimming out Lagrange multipliers is 



H D = J d^[N a je a + NJe) 



(4) 



where N a {a = 1,2,3) is the shift vector, N is the lapse 
function, J€ a is the super-momentum constraint and J€ 
is the super- Hamiltonian constraint. The super-moment- 
um, corresponding to the Oa components of Einstein's 
equations, encodes invariance under spacetime diffeo- 
morphisms within the three-dimensional spatial surfaces 
on which one integrates. The super-Hamiltonian (the 
00 component of Einstein's equations) both encodes in- 
variance under time reparametrizations and generates 
the dynamics (time evolution) of the system. Symmetry 
and dynamics are thus entangled. Canonical quantiza- 
tion follows by promoting the first-class constraints J€ a 
and J€ to operators acting on a Hilbert space of wave- 
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functional yp. Quantum dynamics is then fully speci- 
fied by the equations JfcaV = and the Wheeler-DeWitt 
equation 



(5) 



In a fully background-independent theory, both J€ a 
and jfc are written in terms of the canonical variables as- 
sociated with the fundamental degrees of freedom (met- 
ric and matter) of the system. These expressions are non- 
linear and, in practice, it is extremely difficult to solve the 
constraint equations and construct the physical Hilbert 
space. Symmetry reduction (at the classical level) to the 
flat, homogeneous and isotropic Friedmann-Lemaitre- 
Robertson-Walker (FLRW) metric g MV = {-\,a 2 {t),a 2 {t), 
a 2 {t)) greatly simplifies the problem. The momentum 
constraint is composed only of spatial derivatives and 
it vanishes identically. After integrating over the spatial 
volume (formally divergent but regularizable), the super- 
Hamiltonian in the presence of a matter scalar field ty 
with potential V reads 



1 

2^3 



6k 2 



+ a 



3 K 



(6) 



where = -6aa/ N and = arty I N are the momenta 
conjugate to a and ty, respecively k 2 = 8nG, and K = 0, ± 1 
is the curvature of spatial slices. The constraint 36 = is 
nothing but the first Friedmann equation 



H 2 = — 

3 



ty 2 

Y + v(ty) 



(7) 



The other classical equation of motion is that for the 
scalar field, 



ty + 3Hty+V ( p{ty) = 0. 



(8) 



Quantizing expression (6) and promoting a and ty 
to multiplicative operators and the momenta to deriva- 
tive operators := —i 
J&y[JY,ty] = 0, where 



) a and n := 



id^, one obtains 



„-3.yK 



6 dJY 2 dty 



— +2e GJr V{ty)-%e 



(9) 



and jV = In a is the number of e-foldings. This equation 
may not necessarily be regarded as fundamental. 1 How- 
ever, it gives the correct result in the semi-classical limit, 



Apart from the issue of symmetry reduction, the actual quan- 
tization of the putative full theory can lead to an altogether 



and one can assume it as an effective description of the 
quantum universe in this regime. 

During inflation, the scalar field varies very slowly 
and its kinetic term is negligible with respect to the po- 
tential (slow- roll regime); at the quantum level, it corre- 
sponds to dropping the d 2 ^ term in Eq. (9). Assuming a 

quadratic potential V{ty) = ^m 2 ty 2 , from the Friedmann 
equation (7) it follows that 



6H 2 



m 2 ty 2 . 



(10) 



Thus, the energy scale of inflation sets what shall later 
play the role of quantum correction. 



2.1.2 Perturbations 

When inhomogeneities are switched on, the FLRW mini- 
superspace framework breaks down and one should con- 
sider the full Dirac Hamiltonian (4). Since the super- 
momentum and super-Hamiltonian constraints are non- 
linear in the canonical variables, the problem quickly be- 
comes intractable unless one resorts to some approxima- 
tions. Inflationary inhomogeneous fluctuations are very 
small, so linear perturbation theory is sufficient to ob- 
tain the spectra. The matter scalar is decomposed into a 
homogeneous background (representing the vacuum ex- 
pectation value of the field) and a fluctuation, ty{t,x) = 
</>(£) + Sty[t,x). In this section we ignore the metric back- 
reaction 8g^ v , in which case the scalar is regarded as a 
"test" field. In the standard cosmological model, backre- 
action does not affect the power spectrum at lowest or- 
der in perturbation theory and in the slow-roll trunca- 
tion. This suffices for our purposes also in WDW quan- 
tum cosmology. (We shall include backreaction in the 
LQC case.) The scalar perturbation is decomposed into 
Fourier modes, 



8ty{t,x) = ^Sty k {t)e 
k 



kx 



(11) 



where we assumed spatial slices to be compact (K = 1) 
and the Fourier mode depends only on the modulus k = 



different expression; LQG is an example. Also, canonical quan- 
tum gravity may be embedded in a more general field-theory 
approach such as group field theory [1]. The wave-function <P 
is promoted to a field and the Wheeler-DeWitt equation (5) 
receives non-linear corrections. Linearizing, one gets an effec- 
tive Hamiltonian which can be considered also at the level of 
mini-superspace [10]. 
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|k|. Replacing <p{t) with c/>U,x) in the WDW equation (9), 
the mini-superspace is extended to include also the in- 
finity of modes 8(p k . The wave-function ^[JY,^, {8(p k } k } 
can be actually factorized as a background part times the 
rest, V[Jf, (/>, {6<t> k } k ] = V [^, 4>] Uk>o </>. 8<l> k ] . In 
doing so, one drops self-interaction terms which are 
consistently negligible in first-order perturbation theory. 
Eventually, one obtains [11, 12] 



is Eq. ( 1 6) in the long wave -length limit fc « fc* , then eval- 
uated at k = fc» . This is the n-th order power spectrum 



^ n \k):=^P^\k«K)\ 



(17) 



27r 2> ^—*'\k=k,- 

The lowest-order result coincides with the standard one, 



6 dJ{ 2 



d8<p\ 



e 6 ^m 2 + e^t 



+ e 



8(p 



6H 2 



2 e \2n 



where 



H 

y/ k [^,6(f> k }*0, (12) e:=-— = — 



K 2 <p 2 

TIT 2 



(18) 



(19) 



where y k [Jf,8(p k ] = T'oU/K ,<£]¥ ' k [jV ,(p,8(p k ] and the <p 
dependence is omitted because we used the slow-roll ap- 
proximation (10) to express the background potential in 
terms of the Hubble parameter. 

Noting that JV and 8(j) k correspond, respectively, to 
slow- and fast-evolving variables, at this point one can 
make a Born-Oppenheimer approximation on the solu- 
tion [2, 13]. The latter is written as 



y/ k [Jf,8(p k ] = exp[iS{J/,8(t) k )] 



(13) 



and the functional S is expanded in m 2 x = 3/{4ji£ 2 ] ) = 
6/k 2 : S = m^So + Si + m p 2 S2 + .... Plugging the Ansatz 
(13) into Eq. (12) and expanding, the Oim^) and 0{m 2 ^ 
terms imply So = + e 3J{ H/3, while at the next two orders 
one finds two equations for the wave-functions 



y k l) [jr,8<Pk\ 



(14) 
(15) 



where A and B are chosen to match the amplitudes in the 
WKB approximation. 



2.1.3 Observables 



The wave-functions i//^ 0) and y/^ have been computed 
semi-analytically in [14, 15], to which we refer the reader 
for details. From the explicit solutions, one can calculate 
the two-point correlation function 



— /n,,M 



(16) 



of the scalar perturbation order by order. This quantity 
is directly related to the imprint of inhomogeneous fluc- 
tuations in the cosmic microwave background. However, 
only perturbations which left the comoving Hubble hori- 
zon {aH)~ l =: k~ l and later reentered it can be observed 
in the sky. Therefore, the actual cosmological observable 



is the first slow-roll parameter. Since both H and e are 
approximately constant during inflation, the spectrum 
(which we sometimes call "classical" because of the ab- 
sence of quantum-gravity corrections) is almost scale in- 
variant. 

The next-to-lowest-order expression is the standard 
one times a quantum correction [14]: 



dk)*2?W{k) = 0>W [k]C 2 i 



where 



43.56 H > 



fc 3 



772 



PI/ 



189.18 H 



2 \ 



247.68 H 2 1 

= 1 n T + TZ° 

k 3 ml fc 6 



pi 



fc 3 



in 



(20) 



(21a) 



Pit 



in 



Pl 



l-8 mw {k) + 0{8 : 



2 ) 



(21b) 



and we dubbed the leading Wheeler-DeWitt quantum 
correction 



10 

8 WDW {k):=—{e Pl H) 2 . 



(22) 



C k — ► 1 in the small-scale limit (fc — ► oo), while at large 
scales (fc « fc«) the quantum-corrected power spectrum 
acquires a mild scale dependence which makes the sig- 
nal suppressed with respect to the standard result. A 
similar suppression of the spectrum happens also in 
other models where geometry is quantized, such as non- 
commutative and string inflation [16-19] . At first, it might 
seem counterintuitive that quantum gravity affects large 
scales more than small scales. However, large-scale per- 
turbations left the horizon before (and hence reentered 
after) smaller-scale fluctuations, and they were longer 
exposed to high-energy and high- curvature effects. The 
approximation scheme used to derive Eq. (21a) breaks 
down in the limit C k — » and the critical fc at which that 
happens should not be taken as a physical threshold. 
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From the power spectrum, one can compute the 
scalar spectral index 



1:= 



dln&> & 
dink 



(23) 



which generalizes the definition of an exactly power- 
law type spectrum 3? s ~ k" s ~ l . To calculate this, we no- 
tice that (from aH = k at horizon crossing) dldlnk « 
d/{Hdt) and we recall the background relations, stem- 
ming from the equations of motion, 

2^ 



e = 2He{e 
where 

H(P 



f\ = H{en-K 



77: 



1 (0 







n 



(24) 



(25) 



H 2 [<j>) H 2 (f> 

are the second and third slow-roll parameter, respec 
tively. Since H « const, one gets 

^WDW 



dink 
and 



-35 v 



(26) 



1 = 277-4e + 35 v 



where we have dropped higher-order terms in the com- 
bined <5 WDW /slow-roil expansion. Positivity of the quan- 
tum correction in Eq. (27) ensures suppression of power 
at low wavenumbers. 

The next slow-roll observable is the running of the 
spectral index: 

dn s 



dink 

Combined with Eqs. (26) and (27), it leads to 
a s a 2 [5eri - 4e 2 - £ 2 ) - 95 W dw • 



(28) 



(29) 



The scalar power spectrum expanded to all orders in 
the perturbation wavenumber about a pivot scale fcn is 

ln^ s (fc) = ln^ s {k ) + [n s {k )-l]x+^^-x 2 



no 

m=3 



a 



On) 



(fco) 



nr. 



where x := ln(fc/fc ). As the order of the observables 



„(m) ._ 



d m - 2 a s 
{dink)" 1 - 2 



0{e m )-{-3) m 8 v 



(30) 



(31) 



increases, the classical part becomes smaller and smaller 
but the leading-order quantum correction survives. At 
some order m, the quantum correction will dominate 
over the standard part. Taking (31) into account, Eq. (30) 
can be recast as 



ln0» s (fc) « ln^f (fc) + 5 WDW (fc ) 



(32) 



2.2 Experimental bounds 

Equation (1) is written in units where k is dimensionless. 
In fact, one should make the replacement k — <■ fc/fc m in, 
where fc m i n ~ 1.4 x 10~ 4 Mpc _1 is the largest observable 
scale. Here we used the fact that comoving wavenum- 
bers and multipoles are approximately related by k a 
£/r , where t » 14.4Gpc is the comoving particle hori- 
zon today, and that the lowest early-universe contribu- 
tion to the CMB spectrum is the quadrupole £ = 2. One 
can reexpress 5 WDW in terms of spherical multipoles, and 
k/k m i n = £/£ m i n = (12. A more generous estimate for 
the quantum correction will stem by replacing k m m by 
the pivot scale fcn » fc m in. which we adopt from now on: 
fc — ► fc/fcn. 

The WMAP7 mean for the scalar amplitude in the ab- 
sence of tensor signal is 5^ (fco) = (2.43 ± 0.11) x 10~ 9 at 
fcn = 0.002 Mpc 1 and 68% confidence level (C.L.) [20], 
where the pivot scale fc corresponds to a CMB multi- 
pole £ « 29. Equation (18) and the inflationary condi- 
tion e < 1 yield the upper bound 
-5 



w r/piH) < 9 x 10" 



(33) 



as anticipated in Eq. (1). The bound can be recast for 
the Hubble parameter alone, H < 3.2 x 10 15 GeV, or, 
via the classical equation of motion (7) in the slow-roll 
approximation, for the inflaton potential, V 114 < 6.8 x 
10 16 GeV. In particular, the WDW quantum correction is 
constrained to be 



,(fc ) < 7.9 x 10 



-6 



(34) 



With 

^min instead of fcn the quantum correction is further 
suppressed, 5 WDW (fc ) < 2.6 x 10~ 9 . 

Even taking the upper bound £?\H = 10~ 4 , 5 WDW = 
10~ 5 , quantum corrections are too small to be detected. 
Their dependence on the inflationary energy scale is cru- 
cial for this result. Another reason is that at large scales 
cosmic variance is the leading source of error. The latter 
is a manifestation of the failure of the ergodic theorem 
for the discrete CMB multipole spectrum. For the power 
spectrum {£), cosmic variance is given by [21,22] 

V^ s {£)=^- [ ^ 2 {£). (35) 

Quantum- gravity corrections should be compared with 
the error bars due to cosmic variance with respect to the 
classical spectrum g?®\£). The latter, Eq. (30), is deter- 
mined up to the normalization 8? s {(q), so that the region 
in the (f,^ s (fl/^ s (^)) plane affected by cosmic vari- 
ance is roughly delimited by the two curves 



Var^,o,W 



,(0) 



(£0) 



1 + 



2£+l 



,(0) 



«(0) 



(36) 
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where we take the classical spectrum as reference. The 
WDW-corrected spectrum is given by Eq. (32). In the ab- 
sence of tensor modes and running, the WMAP+BAO+Hn 
dataset (combination of WMAP7 data and observations 
of baryon acoustic oscillations and the Hubble expan- 
sion) yields a scalar spectral index n s (k ) = 0.963 + 0.012 
at k = 0.002 Mpc 1 and 68% C.L. [23]. This number can 
change depending on the priors, but not much. Classi- 
cally, this corresponds to slow-roll parameters at most of 
order e ~ O(10~ 2 ). Therefore, the standard spectrum in 
Eq. (32) can be approximated by \n^°\k) « ln^ s (fco) + 
[« s (fco) ~~ 1]*+ ja s [ko)x 2 . To plot the WDW spectrum, we 
only need to plug in values for the scalar index and its 
running. First, we recast the observables in terms of a set 
of slow-roll parameters dependent on the field potential 
(e.g., [5]): 



£v ■ = 



1 

2k 2 



1 v, 



— I , Vv= 



,00 



K V 



e2 .. 
>V •' 



V,<pV,<p<p<p 

K 4 V 2 ' 



The scalar index (27) and its running (29) become 
n s -l = -6e^ + 2?7 1 /-i-3r5wDw> 

a s = -24e 2 v + \6e v T] v -2^ 2 v -96 WDW . 
For a quadratic potential V(_(f>) oc <p 2 , 
2 



K z (p 



0. 



(37) 

(38) 
(39) 

(40) 



This allows one to reduce the slow-roll parameters to just 
one. A realistic theoretical value for e v at the pivot scale 
isey(fco) = 0.009. 

As shown in Fig. 1, WDW quantum corrections are ex- 
tremely small even in the most generous estimate, and 
they are completely drowned by cosmic variance. 



3 Loop quantum cosmology and 
observations 

3.1 The model 

3.1.1 Homogeneous background 

Loop quantum gravity is based upon a first-order for- 
mulation of gravitational degrees of freedom in terms of 
the spatial densitized triad field E" and the Ashtekar- 
Barbero connection A l a , where i = 1,2,3 is an internal 
index in the su{2) algebra. While the connection is not 
quantized into a well-defined operator, its SU (2) -valued 




500 1000 



Figure 1 (online color at: www.ann-phys.org) Log-linear plot 
of the Wheeler-DeWitt primordial scalar spectrum for 
a quadratic inflaton potential, with e v (ko) = 0.009 and for 
the pivot wavenumber k = 0.002 Mpc -1 , corresponding to 
£ = 29. The shaded region, delimited by the two curves (36), 
is affected by cosmic variance. The inset shows the negligible 
difference between the standard "classical" spectrum (dashed 
line) and the spectrum with Wheeler-DeWitt quantum correc- 
tions (solid line), at 2 < £ < 2.5. 



holonomy along an edge (with representation defined 
by the edge spin label j e ) is a sensible operator. Thus, 
the basic quantities to be quantized are fluxes (integrals 
of the triad on spatial surfaces) and holonomies. In a 
(quasi-)FLRW universe, the densitized triad and connec- 
tion both reduce to one non-trivial component, £ a = 
pSf and A l a = cS' a , where p = a 2 , c = jd, and elementary 
edge lengths are all equal to some common value L = a£o- 
Then, fluxes reduce to F = £^p = L 2 and holonomies 
along an edge of comoving length £q are h = exp(^nT;c) = 
cos(^oc/2)+2t; sin [£oc/2], where t ,■ = io ill are Pauli ma- 
trices. 

We point out that in an exactly FLRW background the 
universe is perfectly homogeneous and there is no mean- 
ingful way to subdivide it into small cells of proper size L. 
Thus, the comoving scale £q is actually arbitrary and cor- 
responds to the size 7^ 1/3 of the fiducial volume in which 
the Hamiltonian constraint is defined. On the other hand, 
in the presence of inhomogeneities the lattice picture 
makes sense (because sub-volumes can be distinguished 
from one another) and fluxes and holonomies can be 
defined on each individual cell, not on the overall fidu- 
cial volume. Quantum dynamics is then expected to have 
the cells vary with time. The freedom to choose a global 
clock in a quasi-homogeneous scenario allows us to pick, 
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e.g., the scale factor a as the time variable, and to regard 
L = L{a) as time dependent. This is the so-called lattice- 
refinement picture [8, 24, 25] . 

A crucial consequence of lattice refinement is that 
inverse-volume corrections are physically meaningful. 
These quantum corrections arise due to the presence 
of inverse -volume expressions in the super-Hamiltonian 
constraint, both in the gravity and matter sector. Inverse 
volumes (i.e., inverse powers of the determinant of the 
densitized triad) are not densely defined operators and 
they must be rexpressed by the so-called "Thiemann's 
trick" in terms of holonomies and positive volume pow- 
ers, at the classical level before quantizing. 

Therefore, contrary to the WDW model, also the back- 
ground equations of motion (and the slow-roll parame- 
ters as well) are deformed by quantum corrections. For a 
matter scalar field, one has 



W = — a 



+ 1 



d\nv 
d\np 



<p + vV <p = Q, 



(41a) 



(41b) 



where a{a) and v{a) are inverse -volume corrections in 
the gravity and matter sector, respectively. Later we shall 
be interested in the semi-classical limit where quantum 
corrections are small, in which case 



a{a) = 1 + a S [nv {a) , v{a) = 1 + v S [nv {a) , 



(42) 



where an and vn are positive constants (calculable in a 
pure FLRW case, arbitrary in the lattice refinement pic- 
ture) and 



fin 



L{a) 



oc a 



(43) 



Here m is an 0(1) constant dependent on the quanti- 
zation scheme (e.g., [25, 26]) and a is determined by a 
power-law Ansatz for the function L{a). While a natural 
value is ex = 6 in pure FLRW, in lattice refinement the only 
constraint is a > 0. In general, however, the background 
inflates ifo-< 0(1) [25]. 

Another type of quantum effect in the dynamics, 
holonomy corrections, is realized in a highly non-linear 
fashion (by construction, from the exponentiation h of 
curvature components) and it becomes important when 
the Hubble radius is about the size of the lattice scale, 
H 1 ~ L. From the classical Friedmann equation H 2 = 
8jiGp/3, this regime heuristically defines the energy den- 
sity (2) and the holonomy correction S^ i := p!p QG . The 
Friedmann equation is modified accordingly but we do 
not show its explicit form here. Further details can be 
found in [5,7,8]. 



3.1.2 Perturbations 

To obtain the dynamics of inhomogeneities, we follow 
the effective constraints method (see [27-32] for other ap- 
proaches). The strategy of applying perturbation theory 
in the classical constraints differs from the one employed 
in standard cosmology (where the action or the Einstein 
equations are perturbed), although it was considered in 
the past [33] . Perturbing the Ashtekar-Barbero variables, 
E" = pSf+SEf, A l a = c8 l a +8A' a , and imposing commuta- 
tion relations among the perturbation components, one 
works out the perturbed form of the seven first-class con- 
straints: the super-Hamiltonian, the three components 
of the diffeomorphism constraint, and the three com- 
ponents of the Gaufi constraint generating infinitesimal 
su{2) gauge transformations in the internal space. To 
capture loop quantum gravity effects, however, one con- 
siders effective constraints encoding inverse -volume 
and/or holonomy corrections. For instance, the effective 
Hamiltonian constraint with inverse-volume corrections 
is assumed to be 



<e[N] 



■ J d 3 xN[a{E)J? g + v{E)J€ n + p(E)J?v + J€v\ , 



(44) 



where N is the lapse function, J€ & , J€ n , JC^j and J€y 
are the contributions of, respectively, gravity, the scalar 
field momentum, spatial Laplacian and potential, and a, 
v and p are correction functions (which depend only on 
the densitized triad [26]). These functions can be taken 
to be of the form 1 + 0(<5i nv ) in the semi- classical limit. 

Closure of the effective constraint algebra must be 
imposed for consistency, {^a,^} = f a ^{A,E)^€ c . The ab- 
sence of anomalies is guaranteed by introducing coun- 
terterms in the algebra (and, hence, in the perturbed 
equations of motion). After some early works based on 
toy models where the constraint algebra was not closed 
explicitly [34-40], the full set of constraints with small 
inverse-volume corrections was derived for vector [41], 
tensor [42], and scalar modes [43,44]. The gravitational 
wave spectra have been studied in [45, 46], while the 
scalar spectrum and the full set of linear-order cosmo- 
logical observables were found in [25] . The observability 
of and experimental constraints on the quantum correc- 
tions were finally considered in [26, 47, 48] . 

The closure of the constraint algebra has been veri- 
fied also in the presence (only) of holonomy corrections 
for vector and tensor modes [41, 42, 49], as well as in the 
scalar sector [50-52]. Cosmological observational signa- 
tures of holonomy effects have been studied for the ten- 
sor sector alone [45, 53-56] . There are indications that 
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holonomy corrections are not significant in the energy 
regime of inflation, but only at near-Planckian densities 
[57]. This is suggested by effective equations for certain 
matter contents with a dominating kinetic energy [58, 59] . 
Another argument is the following [26]. Inverse -volume 
and holonomy corrections are related to each other by 



l 8nG n A~ (pq g ^ 



Pin 



(45) 



Inverse -volume corrections are sizable when the quan- 
tum gravity density (not the inflationary one) is close 
to the Planck density. They can be still large at small 
energy densities, where however holonomy corrections 
are small. Thus, as the energy density decreases in an 
expanding universe there is a competition of the rela- 
tive size of inverse-volume and holonomy corrections, 
the latter falling to small values when the former can 
be still large. For instance, in the inflationary regime (1) 
and for the typical value m = 4 Eq. (45) yields S^i ~ 
1(T 8 / \j <5i nv , and having small holonomy corrections of 
size <5hoi < 10 6 would require inverse -volume correc- 
tions larger than 8 mv > 10~ 4 . 

In the presence of small inverse -volume corrections, 
after anomaly cancellation the system of perturbed equa- 
tions for scalar and tensor modes (vector modes are 
damped during inflation) reduces exactly to two equa- 
tions: 



-A+ — 

z 



u = 0, 



w" + \a 2 A- — \w = 0, 



(46a) 



(46b) 



where primes denote derivatives with respect to confor- 
mal time (' = d T = ad t ), u = z£& is the Mukhanov-Sasaki 
variable encoding scalar perturbations, z = [({>'/ H)[l + 
(an/2 - v n )£i nv ] is a background function (quantum cor- 
rected as well), Si is the gauge-invariant comoving cur- 
vature perturbation (its LQC expression can be found in 
[25,44]), s 2 = 1 + x(«o. vo,CT)5i n v is the square propaga- 
tion speed of the perturbation (discussed in [25]), % = 
(fjvo/3)(o76+ 1) + (an/2)(5 - a/3) (positive in all reason- 
able scenarios), A is the spatial Laplacian, w = ah is 
the gauge -invariant variable associated with both tensor 
modes, and a = a{l - anOi nv /2). 

The parameter space is extended to include the co- 
efficients appearing in Eqs. (42) and (43). However, self- 
consistency of the constraint algebra imposes a condi- 
tion among an, vn and a, thus making one of them de- 
pendent [25]: 



H- v °(H(H =0 - 



The fact that scalar perturbations reduce to just one de- 
gree of freedom u obeying a closed equation is related 
to conservation of 52 at large scales [25] . Failure of clos- 
ing the algebra exactly would immediately spoil also this 
properly. 

The somewhat unexpected possibility that LQC quan- 
tum corrections be large even during inflation is a re- 
flection of the way these corrections enter the physics: 
The structure of spacetime itself is deformed by quan- 
tum effects, via the effective constraints. The theory is 
diffeomorphism invariant, but not with respect to the 
standard classical transformations. Gauge transforma- 
tions belonging to a deformed algebra no longer corre- 
spond to ordinary coordinate transformations on a man- 
ifold. Thus, to take the new gauge structure into account 
one should rely only on gauge-invariant perturbations. 
This philosophy (first quantize the classical system, then 
cast it in gauge -invariant variables) is embodied in the 
Mukhanov equations (46). 

One might wonder whether one would get the same 
results by fixing the gauge before quantizing. However, 
gauge fixing and quantization do not commute because 
the latter deeply affects the very notion of gauge in- 
variance. Whenever gauge-ready variables can be con- 
structed after quantizing, the gauge -invariant approach 
must be preferred. This is well known also in standard 
cosmology [60] . The price to pay in doing otherwise is, 
in the worst case, to produce unphysical perturbative 
modes, or at best to obtain an incomplete version of 
the perturbed quantum equations. Also ignoring back- 
reaction of the metric and considering just a perturbed 
test scalar is undesirable, contrary to the WDW case, 
because backreaction contributes to the actual form of 
quantum gauge transformations and hence of the gauge- 
invariant variables. Again, this can lead to an incomplete 
treatment in partial disagreement with the full gauge- 
invariant equations. 



3.1.3 Observables 

The scalar spectrum is the expectation value of S/i over a 
momentum ensemble at large scales, evaluated at hori- 
zon crossing: 

fc 3 



(\u k «kj > 



k=k, 



(48) 



Solving the Mukhanov equation (46a) asymptotically and 
plugging the solution in the above formula, one obtains 
the LQC version of Eq. (20) with 

o-a % 



(47) q«i + rA 



7s:=v (| + l J + . 2( 



(7 + 1 



(49) 
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In the limit case a = 0, the quantum correction is con- 
stant and the only change with respect to the classical 
case is the normalization of the spectrum. Then, y s = 
vn - 5ao/2 could be of either sign. If a ^ 0, there is a 
large-scale enhancement of power because 8 mv ~ k~ a at 
horizon crossing and y s > due to the dominating term 
oc r 1 . Similarly, the scalar index is 

n s -1^2r]-4e+aj ns S inv , y„ s := a -2v + -, (50) 

<7+ 1 



while the scalar running reads 



a s a 2(5er? - 4e 2 - £ 2 ) + cr(4e - aj ns )S inv , 



(51) 



where e := a {al2 + 2e-rfj + v (cr/6- l)e. 

Due to the possibly large size of the quantum correc- 
tions, it will be useful to complete the set of first-order 
observables and include also the tensor sector. The grav- 
itational spectrum is 



32G k 3 , n, . 

leading to [25,46] 
^ t = 64^G|£) (1 + yAnv). 



y t := 



(T-l 
CT+1 



a - 



The tensor index is 
din &> x 



n t :-- 



2e-o-y,5 il 



d\nk 

Finally, the tensor-to-scalar ratio is 

r := — « 16e[l + (y, - y s )«5 inv ] , 

^ s 

which yields the consistency relation 
r = -8{n, + [n t (y t - y s ) + <ry t ]<W , 



(52) 



(53) 



(54) 



(55) 



(56) 



to be plugged into numerical codes in the place of the 
classical one r = -8n t . 



3.2 Experimental bounds 

Because of the delicate interplay between quantum cor- 
rections and the requirement of intersecting the allowed 
windows in the parameter space in a common consistent 
region, the possibility clearly arises that this loop model 
of quantum cosmology be falsifiable by near-future ob- 
servations. The present status at least provides stringent 
bounds on quantum corrections. 

As in section 2.2, one rewrites the observables in 
terms of the potential-dependent slow-roll parameters 



(37); the resulting lengthy expressions can be found in 
[26]. Since 



m„m-\ 



/s<5inv(fco), 



where 



/.:= 



o-[3a (13o--3) + v (r(6+ 11a)] 

18(CT+1) 



(57) 



(58) 



the scalar spectrum (30) becomes 
9» a m ~ ^s(fco)exp|[« s (A:o) 



l]x+ X 



+ fs8invik ) 



x\l 



-ax 



1 

+ -{e~ 
a 



1) 



(59) 



This is the expression to be used in numerical analyses 
and when comparing the LQC signal with cosmic vari- 
ance. 

Before doing so, we notice the existence of a theo- 
retical upper bound on the quantum correction 5 LQC := 
ao^inv (Equation (47) allows to remove vn from parame- 
ter space, except in the case a = 3 which can be treated 
separately.) For the validity of the linear expansion of the 
perturbation formula? where the 0(5i nv ) truncation has 
been systematically implemented, we require that 



'LQC 



-f) =5 LQ c(fco) y <1 (60) 



for all wavenumbers relevant to the CMB anisotropics. 
For the pivot scale £o = 29, the quadrupole £ = 2 gives 
the bound S LQC {k ) < 8™™ = 14.5 _<J , shown in Tab. 1 for 
some choices of a. 

To illustrate some of the possibilities CMB data ma- 
nipulations can offer to constrain quantum gravity mod- 
els with free parameters, we recall the likelihood analysis 
carried out in [26, 47] for the quadratic potential (among 
others). The Cosmological Monte Carlo (CosmoMC) code 
[61] was run with the data of WMAP7 [23] combined with 
large-scale structure (LSS) [62] (including BAO), HST [63], 
Supernovae type la (SN la) [64] , and Big Bang Nucleosyn- 
thesis (BBN) [65], assuming a ACDM model. Figure 2 
shows an example of likelihood profile for a = 3/2 in 
the plane (e v ,8 LQ c)- Both parameters are evaluated at 
the pivot scale fcn = 0.002 Mpc -1 . Obviously, negligible 
or exactly vanishing quantum corrections are compati- 
ble with observations. On the other hand, from the 95% 
confidence-level contour one sees that quantum correc- 
tions above 8 LQC [ko) > 1.7 x 10~ 3 can be excluded. This 
and the upper bounds for various a's are reported in Tab. 
1. Except for extreme values a « 1, the observational up- 
per bounds are consistent with the theoretical prior, thus 
verifying an important internal check of the model. 
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Table 1 Theoretical priors on the upper bound <5™^f and 95% C.L. upper 
limits of S LQC = ao<5i nv constrained by observations for a quadratic poten- 
tial with different values of a and at the pivot scale k = 0.002 Mpc -1 [ ] 
The likelihood analysis is omitted for a = 6 since the signal is below the 
cosmic variance threshold already when a = 2. For a = 3, the parameter 
8 lqc = vnSinv has been used. 



a 


0.5 


1 


1.5 


2 


3 


6 


jrmax 

U LQC 


0.26 


6.9 x 10" 2 


1.8 x 10" 2 


4.7 x 10" 3 


3.2 x 10" 4 


1.0 x 10" 7 




0.27 


3.5 x 10 -2 


1.7 x 10 -3 


6.8 x 10" 5 


4.3 x 10 -7 





. x 1 0"' 



1.5 - 




0.015 



0.02 



Figure 2 (online color at: www.ann-phys.org) Two- 
dimensional marginalized distribution for the inverse-volume 
LQC quantum correction <S LQC {k ) and the slow-roll parame- 
ter e v {ko) with the pivot fco = 0.002 Mpc -1 for a = 1.5 and 
a quadratic potential, constrained by the joint data analysis 
of WMAP7, LSS (including BAO), HST, SN la, and BBN. The 
internal and external lines correspond to the 68% and 95% 
confidence levels, respectively [26]. 




500 1000 



Figure 3 (online color at: www.ann-phys.org) Log-linear plot 
of the LQC primordial scalar spectrum g? s {_£) with inverse- 
volume quantum corrections for a quadratic inflaton potential, 
with e v (Ar ) = 0.009 and for the pivot wavenumber k = 0.002 
Mpc -1 , corresponding to ( = 29. The classical case is rep- 
resented by the dotted line, while solid curves correspond to 
a = 1,1.5,2 (decreasing thickness). The shaded region is af- 
fected by cosmic variance. 



4 Non-Gaussianity 



Comparing the table entries with the upper bound for 
the WDW quantum correction, Eq. (34), we see that LQC 
inverse -volume corrections can be orders of magnitude 
larger when a < 2. The scalar power spectrum for vari- 
ous values of a is shown in Fig. 3 against cosmic variance. 
When a < 1, quantum corrections are strong enough to 
overcome the error from cosmic variance. Whether these 
parameter values are realistic in a more complete theory 
remains, however, to be seen. 



The effect of quantum corrections goes beyond linear 
perturbation theory and higher-order observables can 
be calculated. As the perturbative level is increases, the 
statistics of inhomogeneous fluctuations deviates from 
the Gaussian one and odd-order correlation functions ac- 
quire non-vanishing values. In particular, the bispectrum 
(three-point correlation function of the curvature pertur- 
bation) can be constrained by observations. 

To the best of our knowledge, there is only one work 
on inflationary non-Gaussianiry in loop quantum cos- 
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mology [66], and none in the WDW case. A detailed 
calculation of the momentum-dependent bispectrum 
shows that no appreciable LQC signal can be detected. 
We can in fact reach the same conclusion here by a 
model-independent shortcut, valid only in the so-called 
squeezed limit (constant non-linear parameter) but be- 
yond perturbation theory and both for LQC and WDW 
quantum cosmology. 

Let ( be the curvature perturbation on uniform den- 
sity hypersurfaces. The latter is a gauge-invariant quan- 
tity proportional to the comoving curvature perturbation 
5? in standard inflation; their relation in the presence of 
inverse-volume corrections has not been studied yet, but 
what follows is fairly independent on this detail. In mo- 
mentum space, the three-point correlation function of ( 
is 

<fk/k 2 Ck 3 > =: (27r) 3 <5(k! + k 2 + k. 3 )B ( (k lt k 2 , k 3 ) , (61) 
where B^ is called bispectrum and it is defined by 

B f (fci.fcz.fca) = \fm.Vcx,k 2 ,h) £ P ( {k a )P ( {kp) , (62) 

where a,/3 = 1,2,3, /nl is called non-linear parameter 
and is momentum dependent in general, and is the 
spectrum of ( . The form of the non-linear parameter de- 
pends on the model of primordial perturbations. In the 
simplest case [67-69], one decomposes the non-linear 
curvature perturbation £nl(x) into a Gaussian linear part 
( and a non-linear part: 



C NL = C + C N = C+^° L cal (C 2 -(C 2 » ) 



(63) 



where the non-linear parameter f^ al is constant. By def- 
inition, (Cnl) = (0 = 0. Then, a direct calculation of the 
bispectrum shows that 



fNLiki,k 2 ,k 3 ) = f^ ai . 



(64) 



In fact, the Fourier transform of the non-linear part is 



local 



>-i\ _ _ f J 
''k ~ 5 -/NL 



-(27r) J <5(k)<0 + 



/ 



d 3 p 



CpCp k 



(65) 



The first term stems from the fact that the auto -correlation 
function is x independent. Since all momenta must not 
vanish at the same time, this piece can be thrown away. 
The second term enters into the three-point function, 
which at lowest order is (e.g., [5]) 



<Ck 1 Ck 2 c k ;> + (k3-k 2 ) + (k3-k 1 ) 



(2^) 3 <5(k 1 + k 2 - kg) -$£*2P ( (fci)P f (fc 2 ) 
5 



+ (k 3 -k 2 ) + (k 3 -k 1 ), 



(66) 



which yields Eq. (64) after comparing Eqs. (61) and (62). 
The decomposition (63) is pointwise in configuration 
space and for this reason it is called local model. For a 
power-law scalar spectrum 3? s oc fc" s ~\ the local bispec- 
trum reads 



B l ° cal {k l ,k 2 ,k 3 ) = -f^A 2 ( Y l 

5 a<p 



{k a kp) 



4-n s 



(67) 



where A( is a constant amplitude. This expression can be 
converted into one with spherical multipoles. 

The expression (67) peaks at the squeezed limit where 
one of the edges of the triangle (ki,k 2 ,k3) collapses [70, 
71]: 



fci a fc? » k 3 



k~0. 



(68) 



Sending, e.g., k 3 — <■ 0, by conservation of momenta one 
has ki — k 2 and 



B l ° ca[ (h, h, k 3 - 0) = y /i° L cal P f (h)P ( (fc 3 ) . (69) 



/k-x 



Measuring the bispectrum in this configuration, one can 
obtain an estimate of £[S . In the local bispectrum, 
small- and large-scale modes are coupled together. 

The squeezed limit can be understood in a fairly intu- 
itive way in all models where the curvature perturbation 
( is constant at large scales [70, 72] . Split ( into a corse- 
grained and a fine-grained perturbation, 

r d 3 k ;kv f d 3 k 

Jk<k, {2ny Jk>k t {2n) 6 

=:Cc(T,x)+C q (T,x). (70) 

In the limit (68), is larger than the Hubble horizon and 
can be treated as constant in time. Then ( (x 3 ) ~ £c(x 3 ) de- 
fines a new coordinate background x' « [1 + £ c (x 3 )]x in- 
side the horizon. In the new coordinates and up to linear 
order, 

C q (x') a C q (x) + (x'-x)~( q (x) 



d 

~ C q (x) + Cc(X3)x- ^~Cq(x) • 



(71) 



If the linear perturbation f q (x) is Gaussian, in the squeezed 
limit we have 

(C( Xl )C(X 2 )C(X3)> ~ (CqftiKqfriKcCxa)) 



« (Cc(X 3 )X 1 ~ [C q (X!)C q (X 2 )] 



= (C 2 (X3)) cXl - — <C q ( Xl Kq(X 2 )) q 



c(0, 



(72) 
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where in the second line we exploited translation invari- 
ance, in the last line we used g = |xi -X2I and dg/dxi = 
xi Iq, and we denoted with 4 the two-point correlation 
functions of (. Since the latter goes as ffp oc p _( " s_1) , one 
gets 



(ffciKteKteO) « -in s - 1)4° (0)4° (£0 - 



(73) 



Comparing this expression with Eq. (69), we finally ob- 
tain 



/■local _^ fi 

/nl ~12 U " 



n s ). 



(74) 



For spectra which are almost scale -invariant [n s - 1 
small) at large scales, the level of non-Gaussianity is very 
low, /nl « 1- Tensor modes produce an even lower signal. 
This result [70, 72] is general enough to be applied both to 
WDW and loop quantum cosmology, which we have seen 
to be compatible with almost scale invariance. Therefore, 
considering the current 95% C.L. bound on the local non- 
linear parameter coming from combined CMB and large- 
scale structure [73] observations, -5 < /^° cal < 59 [23], 
the non-linear parameter in the squeezed limit is small 
and the quantum corrections considered here have no 
appreciable impact on the bispectrum. 



5 Outlook 

Quantum gravitational effects modify the spectra of cos- 
mological perturbations and their imprint in the cos- 
mic microwave background. In this review, we compared 
two canonical approaches, the one based on the usual 
Wheeler-DeWitt quantization and loop quantum cos- 
mology. Wheeler-DeWitt quantum corrections are too 
small to be detected, even in the most optimistic upper 
bound, Eq. (34). The model therefore is not falsifiable, at 
least under the assumptions made in the derivation of 
the results, but at least it is compatible with what we ob- 
serve. 

In contrast, LQC inverse -volume corrections can be of 
much greater size and produce an enhancement, rather 
than suppression, of the large-scale spectra. While in 
the WDW case quantum corrections change the inhomo- 
geneous dynamics but leave homogeneous background 
equations unmodified, in LQC the latter are deformed, 
too. However, this is not the reason why LQC effects are 
potentially several orders of magnitude larger than the 
WDW quantization. Rather, the key ingredient is the scale 
compared with the Planck energy density p P i in the ratio 
defining the quantum correction: for WDW it is the in- 
flationary scale p in fl, for LQC it is determined by the char- 



acteristic discreteness scale of the semi- classical state de- 
scribing the quantum universe. This effective energy den- 
sity can be as large as the Planck density, pi n fl « p QG < 

This also highlights the different origin of the obser- 
vational bounds presented above. While the WDW quan- 
tum correction (34) is constrained somewhat indirectly 
via the usual bounds on the inflationary energy scale, 
in LQC we have some free parameters on which we 
have little control theoretically, due to the formidable 
(and yet unsurmounted) difficulties in explicit construc- 
tions of cosmological semi-classical states in the full 
theory. LQC inverse-volume corrections depend on a 
phenomenological quantum-gravity scale as well as on 
partly heuristic, partly quantitative arguments indicating 
how to implement discrete quantum geometry in a quasi- 
homogeneous cosmological setting. A multi-variate like- 
lihood analysis involving all the cosmological parame- 
ters, including LQC ones, is thus more adequate to the 
task. 

Observations constrain LQC inverse -volume quan- 
tum corrections below their theoretical upper bound, but 
in some instances the signal is above the threshold of cos- 
mic variance. Experiments such as PLANCK or of the next 
generations should then be able to reach the sensitivity 
to detect a quantum gravity signal or, in its absence, place 
yet more stringent constraints. In turn, pressure from ac- 
tual data will stimulate the quest for a better understand- 
ing of the fundamental properties of the states of the full 
theory, and a greater control over parameters which, as 
the discreteness scale L, are presently phenomenologi- 
cal. 
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